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Abstract 

Final results from the DELPHI Collaboration on the lifetime of B"^ and 
B° mesons and the mean b-hadron lifetime, are presented using the data col- 
lected at the Z° peak in 1994 and 1995. Elaborate, inclusive, secondary ver- 
texing methods have been employed to ensure a b-hadron reconstruction with 
good efficiency. To separate samples of B+ and B° mesons, high performance 
neural network techniques are used that achieve very high purity signals. The 
results obtained are: 

Tb+ = 1.624±0.014 (stat)±0.018 (syst) ps 
tbo = 1.531±0.021 (stat)±0.031 (syst) ps 
^ = 1.060±0.021 (stat) ±0.024 syst 

and for the average b-hadron lifetime: 

Tb = 1.570±0.005 (stat)±0.008 (syst) ps. 



(Accepted by Eur. Phys. J.) 
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1 Motivation and Overview 

In addition to testing models of b-hadron decay, knowledge of b-hadron lifetimes is of 
importance in the determination of other Standard Model quantities such as the CKM 
matrix element Kb and in measurements of the time dependence of neutral b-meson 
oscillations. 

The Spectator Model provides the simplest description of b-hadron decay. Here, the 
lifetime depends only on the weak decay of the b-quark with the other light-quark con- 
stituent (s) playing no role in the decay dynamics. This, in turn, leads to the prediction 
that all b-hadron species have the same lifetime. Non-spectator effects however such 
as quark interference, W exchange and weak annihilation can induce lifetime differences 
among the different b-hadron species. Models of b-hadron decay based on expansions 
in 1/mb predict that, T(Ab) < t(B°) ~ t(Bs) < r(B"'"),^ and this lifetime hierarchy has 
already been confirmed by experiment. There is a growing consensus between models 
that a difference in lifetime of order 5% should exist between the B"*" and B° meson [1] 
and it is currently measured to be, ^ = 1.085 ± 0.017 [2]. Clearly, more precise life- 
time measurements of all b-hadron species are valuable in order to test developments in 
b-hadron decay theory. 

This paper reports on the measurement of the B+ and B*^ meson lifetimes from data 
taken with the DELPHI detector at LEP I, in sub-samples separately enriched in B"*" and 
B° mesons. In addition, a measurement of the mean b-hadron lifetime Tb (i.e. with 
Ab, B°, Bg, B"*" unseparated) is obtained, which is a quantity of importance for many b- 
physics analyses at LEP e.g. in the extraction of the CKM matrix element |Kb|- 

The analysis proceeds by reconstructing the proper time of b-hadron candidates (t = 
Lmo/pc where L,p and mo are the decay length, momentum and rest mass respectively) 
and fitting the t distribution in simulation to the data distribution in a ^^-minimisation 
procedure. In the fit for the B"*" and B° lifetimes, the Bs and Ab lifetimes are set to 
their measured values whereas in the fit for Tb, all B-species are assumed to have the 
same lifetime and Tb is the only free parameter in the fit. The approach used is highly 
inclusive and based on the DELPHI inclusive B-physics package, BSAURUS [3] . Aspects 
of BSAURUS directly related to the analysis are presented in a summarised form but the 
reference should be consulted for full details of the package. 

After describing parts of the DELPHI detector essential for this measurement in Sec- 
tion 2, the data sets are described in Section 3 and relevant aspects of BSAURUS are 
highlighted in Section 4. Section 5 describes the reconstruction of the B-candidate proper 
time from measurements of decay length and momentum. Samples with a ~ 70% purity 
in B+ or B° mesons were achieved by use of a sophisticated neural network approach 
which is described in Section 6. Section 7 shows the results of the lifetime fits and finally, 
systematic uncertainties on the measurements are dealt with in Section 8. 

2 The DELPHI Detector 

A complete overview of the DELPHI detector [4] and of its performance [5] have been 
given in detail elsewhere. This analysis depends crucially on precision charged particle 
tracking performed by the Vertex Detector (YD), the Inner Detector, the Time Projection 
Chamber (TPC) and the Outer Detector. A highly uniform magnetic field of 1.23 T 
parallel to the electron beam direction, was provided by the superconducting solenoid 

'^Note that the corresponding charge conjugate state is always implied throughout this paper. B*^ always refers to the 
state and A;, refers to any, weaJfly decaying, b-baryon. 
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throughout the tracking volume. The momentum of charged particles was measured with 
a precision of (Tp/p < 1.5% in the polar angle region 40° < 9 < 140 and for p < 10 GeV/c. 

The VD was of particular importance for the reconstruction of the decay vertices of 
short-lived particles and consisted of three layers of silicon microstrip detectors, called 
the Closer, Inner and Outer layers, at radii of 6.3 cm, 9 cm and 11 cm from the beam line 
respectively. The VD was upgraded [6] in 1994 and 1995 with double-sided microstrip 
detectors in the Closer and Outer layers providing coordinates^ in both rcj) and rz. For 
polar angles of 44° < 9 < 136°, a track crosses all three silicon layers of the VD. The 
measured intrinsic resolution is about 8 /xm for the r0 coordinate while for rz it depends 
on the incident polar angle of the track and reaches about 9 /xm for tracks perpendicular 
to the modules. For tracks with hits in all three r^) VD layers the impact parameter 
resolution was cr^^ = ([61/(p sin^/^^)]^ -|- 20^) //m^ and for tracks with hits in both rz 

layers and with 9 ^ 90°, a^^ = ([67/(p sin^^^ 9)^ + 33^) /im^. Before the start of data 
taking in 1995 the ID was replaced with a similar device but with a larger polar angle 
coverage in preparation for LEP 2 running. The impact of this change on the current 
analysis is relatively minor. 

Calorimeters detected photons and neutral hadrons by the total absorption of their 
energy. The High-density Projection Chamber (HPC) provided electromagnetic calorime- 
try coverage in the polar angle region 46° < 9 < 134° giving a relative precision on the 
measured energy E of as/E = 0.32/Ve © 0.043 {E in GeV). In addition, each HPC 
module is essentially a small TPC which can chart the spatial development of showers 
and so provide an angular resolution better than that of the detector granularity alone. 
For high energy photons the angular resolutions were ±1.7 mrad in the azimuthal angle 
4> and ±1.0 mrad in the polar angle 9. 

The Hadron Calorimeter was installed in the return yoke of the DELPHI solenoid and 
provided a relative precision on the measured energy of as/E — 1.12 /y/E © 0.21 (£■ in 
GeV). 

Powerful particle identification (see Section 4.4) was possible by the combination of 
dE/dx information from the TPC (and to a lesser extent from the VD) with information 
from the DELPHI Ring Imaging CHerenkov counters (RICH) in both the forward and 
barrel regions. The RICH devices utihsed both liquid and gas radiators in order to 
optimise coverage across a wide momentum range: liquid was used for the momentum 
range from 0.7 GeV/c to 8 GeV/c and the gas radiator for the range 2.5 GeV/c to 25 
GeV/c. 

3 Data Samples 

3.1 Hadronic Event Selection 

Hadronic Z'^ decays were selected by the following requirements: 

• at least 5 reconstructed charged particles, 

• the summed energy of charged particles (with momentum > 0.2 GeV/c) had to be 
larger than 12% of the centre-of-mass energy, with at least 3% of it in each of the 
forward and backward hemispheres defined with respect to the beam axis. 

Due to the evolution of the DELPHI tracking detectors with time, and run-specific de- 
tails such as the RICH efficiency, the data were treated throughout this analysis as two 

'^DELPHI has a cylindrical polar coordinate system. The z-axis points along the beam direction and r and if) are the 
radial and aaimuthal coordinates in the transverse plane. 
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independent data sets for the periods 1994 and 1995. The hadronic event cuts selected 
approximately 1.4 million events in 1994 and 0.7 million events in 1995. 

3.2 Event Selection 

Hadronic events were enhanced in Z° — > bb events, by cutting on the DELPHI com- 
bined b-tagging variable described in [7] . In the construction of the b-tag, the following 
four variables were selected that are highly correlated with the presence of a b-hadron 
but only weakly correlated between themselves: 

1. The mass of particles included in a reconstructed b-hadron secondary vertex. 

2. The probability that if a track originated from the primary vertex, it would have 
a positive lifetime-signed impact parameter, with respect to this vertex, at least as 
large as that observed. 

3. The fraction of the total jet energy contained in the tracks associated with a sec- 
ondary vertex, fitted with an algorithm run only on the tracks associated with the 
jet. 

4. The rapidity, with respect to the jet axis direction (see Section 4.1), of tracks included 
in the secondary vertex. 

These variables were combined in likelihood ratios (which assumes they are independent), 
to give a single b-tag variable that can be applied to tag single jets, hemispheres or the 
whole event. The event b-tag, comparing data and simulation, is shown in Figure 1 
together with the contributions of b- and u,d,s,c-events from the simulation. The arrow 
indicates the position of the cut made to enhance the sample in Z° — > bb events. 




Figure 1: The event b-tag variable in 1994 data and simulation, normalised to the number 
of entries. 
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In addition to the b-tag requirement it was demanded that events be well contained 
in the barrel of the DELPHI detector by making the following cut on the event thrust 
axis vector: \cos{9thrust)\ < 0.65. After applying the event selection cuts the purity in 
Z° — > bb events was about ~ 94% and the data samples consisted of 285088(136825) 
event hemispheres in 1994(1995) respectively. 

3.3 Simulated Data 

Simulated data sets were produced with the JETSET 7.3 [8] package with tunings 
optimized to DELPHI data [9] and passed through a full detector simulation [5]. Both 
Z° — > qq and Z° — > bb events were used and separate samples were produced for com- 
parison with 1994 and 1995 data. The same hadronic and event selection criteria were 
applied to the simulation samples as for the data resulting in 640888(239061) events of 
qq Monte Carlo in 1994(1995) and 1581499(422304) events of bb Monte Carlo. 

4 General Tools 

The identification of b-hadron candidates and the reconstruction of their decay length 
and momentum (or equivalently their energy), was made in a completely inclusive way; 
i.e. the analysis was sensitive to all b-hadron decay channels. This section describes 
briefly some tools that were essential in being able to decide whether the data were 
likely to contain a b-hadron and what the probability was, for that b-hadron, to be a 
B+,B'^,Bg or b-baryon. Two specially constructed neural networks (the TrackNet and the 
BD-Net, described below) made it possible to distinguish whether a particle was likely 
to have originated from the primary vertex (in the fragmentation process or the decay 
of an excited state), from the weakly decaying b-hadron secondary vertex or even from 
the B — > D cascade tertiary vertex. The reconstruction of the b-hadron secondary vertex 
position was essential in order to reconstruct the decay length, which in turn tags the 
presence of a long-lived b-hadron state. Exploiting the excellent particle identiflcation 
ability of the DELPHI detector was another aid to b-hadron reconstruction. The presence 
of high-momentum electrons or muons is a powerful signature of a b-hadron decay and the 
detection of kaons or protons also provides valuable information on the type of b-hadron. 

4.1 Rapidity 

Events were split into two hemispheres using the plane perpendicular to the thrust 
axis. A reference axis was defined in each hemisphere along a jet reconstructed via 
the routine LUCLUS [8] with p± as a distance measure and the cutoff parameter 
'^iom(PARU(44))=5.0 GeV/c. In simulation studies this was found to give the best re- 
construction of the initial b-hadron direction. For hemispheres where two or more jets 
were reconstructed (about 16% of the cases), the b-tag applied at the jet level was used 
to discriminate the b-jct from the gluon jet. With this scheme, the probability to select 
correctly the two b-jets in a three-jet event was about 70%. 

The rapidity, with respect to the reference axis, was defined as 

y = i-ln((£; + P||)/(£;-P||)), 

where identified particles were assigned their respective masses and all others were as- 
signed a pion mass. Figure 2, shows the rapidity distribution for real and simulated data 
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after the event selection cuts have been apphed. Particles originating from the decay 
chain of a b-hadron are seen to have higher mean rapidities than particles originating 
from the primary vertex. This property made the rapidity a useful input quantity to the 
TrackNet (Section 4.2) and for b-hadron decay vertex reconstruction (Section 4.3). In 
addition, a Rapidity Algorithm was defined that summed {pi, Ei, mi) for particles with 
7/ > 1.6 in order to form an estimate of the weakly decaying b-hadron four- vector. 



25000 
20000 
15000 
10000 



5000 



°0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Rapidity 

Figure 2: Distributions of rapidity from 1994 data and simulation where the normalisation 
is to the number of entries. From the simulation, tracks originating at the primary vertex 
and tracks from the b-hadron decay chain are also plotted. 



4.2 The TrackNet 

The TrackNet is the output of a neural network which supplied, for each track in the 
hemisphere, the probability that the track originates from the decay of the b-hadron. 
The network relied on the presence of a reconstructed secondary vertex (see Section 4.3) 
providing an estimate of the b-hadron decay position in the event hemisphere. The b- 
hadron four-vector was reconstructed by the Rapidity Algorithm described above. The 
most important inputs to the neural network were: 

• the magnitude of the particle momentum in the laboratory frame, 

• the magnitude of the particle momentum in the B-candidate rest frame, 

• the track hclicity angle, defined as the angle between the track vector in the B- 
candidatc rest frame and the B-candidate momentum vector in the lab frame, 

• a flag to identify whether the track was in the secondary vertex flt or not, 

• the primary vertex track probability (defined as in Point 2. of Section 3.2), 

• the secondary vertex track probabihty (defined as in Point 2. of Section 3.2), 

• the particle rapidity. 

Figure 3, shows the TrackNet distribution for real and simulated data after the event 
selection cuts have been applied. The simulation distribution is divided into contribu- 
tions from Z° — > bb events representing ~ 94% of the total and the remaining 6%, shown 
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in black, is due to tracks in u,d,s-quark and c-quark Z decays. The b-events are fur- 
ther divided into 'signal' and 'background'. The signal consisted of tracks originating 
from the b-hadron weak decay chain and this class of track formed the target for the 
neural network training. The background consisted of everything else in b-events such as 
fragmentation tracks and decay products of excited b-hadrons. Figure 3 illustrates the 
excellent separation, of b-hadron decay products from other tracks, possible using the 
TrackNet and shows a good overall agreement in shape between data and simulation. As 
a result of the b-tagging cut, the contribution from u,d,s,c-events is dominated by c-events 
(~ 80%) which account for the spike around TrackNet= 1.0 in the u,d,s,c distribution. 




TrackNet 

Figure 3: Distributions of the TrackNet from 1994 data and simulation, normalised to 
the number of entries. 



4.3 Secondary Vertex Finding 

For each hemisphere an attempt was made to fit a secondary vertex from charged 
particle tracks that were likely to have originated from the decay chain of a weakly 
decaying b-hadron. Tracks were first selected by imposing the following set of standard 
quality cuts: 

• impact parameter in the rcj) plane \5r-4)\ < 4.0 cm, 

• impact parameter in the z plane l^^l < 6.0 cm, 

• |cos6'| < 0.94, 

• ^ < 1.0, 

• at least one r — (j) track hit from the silicon vertex detector (VD), 
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• tracks must not have been flagged as originating from interactions witli detector ma- 
terial by the standard DELPHI interaction vertex reconstruction package, described 
in [10]. 

In addition, all charged particles were required to have rapidity greater than the value 1.6, 
giving good discrimination between tracks orginating from the fragmentation process and 
those originating from the b-hadron decay chain (see Figure 2). In order to be considered 
further for vertex fitting, charged particles were selected using a procedure based on: 

• The likelihood to be an electron, muon or kaon. 

• The 'lifetime' content of the track based on the three-dimensional crossing point of 
the track with the estimated b-hadron direction vector from the Rapidity Algorithm. 
The direction vector was forced to pass through the event primary vertex position. 

• Particle rapidity. 

If at least two tracks are selected, these tracks are fitted to a common secondary vertex in 
three dimensions with the constraint that the secondary vertex momentum vector must 
point back to the primary vertex error ellipse. If the fit did not converge^ the track 
making the largest contribution was stripped away, in an iterative procedure, and the 
fit repeated. In a final step, tracks with a large TrackNet value, which had not already 
been selected by the initial track search, were added to the vertex track list and the vertex 
re-fitted. 

The secondary vertex described above refers to the 'standard' vertex fit which was an 
essential input to the calculation of other b-physics quantities used in the analysis e.g. 
the BD-Net described in Section 4.5. More sophisticated secondary vertexing algorithms 
developed specifically for the extraction of B lifetimes are described in Section 5.2. 

4.4 Particle Identification 

The MACRIB package [11] provided separate neural networks for the tagging of kaons 
and protons which combined the various sources of particle identification in DELPHI. An 
efficiency for the correct identification of of 90% (70%) was attained with a contamina- 
tion of 15%(30%) forp < 0.7 GeV/c {p > 0.7 GeV/c). The corresponding contamination 
for proton identification at the same efficiencies was, 2%(40%) forp < 0.7 GeV/c (p > 0.7 
GeV/c). 

Electron and muon candidates were defined according to the standard DELPHI Icpton 
identification criteria. Only muon and electron candidates with energy larger than 3 GeV 
were selected. 

4.5 The BD-Net 

In selecting tracks for inclusion in the B secondary vertex fit, there is inevitably some 
background from tracks that originate not from the B decay vertex directly, but from 
the subsequent D cascade decay. When such tracks are present the vertex is in general 
reconstructed somewhere in the region between the B decay point and the D vertex and 
the reconstructed decay length will be biased to larger values. The resolution of the B 
decay vertex is therefore improved if these tracks can be identified and removed from the 
secondary vertex. In order to identify these tracks a neural network (the BD-Net) was 
developed based on the following discriminating variables: 

"^Non-convergence means; the was above 4 standard deviations at the end of the first 10 iterations or above 3 standard 
deviations at the end of the next 10 iterations or took more than 20 iterations in total. 
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• the angle between the track vector and the reconstructed B flight direction, 

• the primary vertex track probability (defined as in Point 2. of Section 3.2), 

• the secondary vertex track probabihty (defined as in Point 2. of Section 3.2), 

• the momentum and angle of the track vector in the B rest frame, 

• the TrackNet output defined in Section 4.2, 

• the kaon network output, described in Section 4.4, 

• the lepton identification tag, mentioned in Section 4.4. 

The network was trained to recognise tracks originating from the decay chain 
B — > D — > X ('signal') compared to all other tracks in b-events ('background') where, 
in addition, all tracks must have TrackNet values larger than 0.5. Figure 4 shows the 
BD-Net variable for data and simulation after the event selection cuts have been applied 
plus the same track selection cuts that are used in the vertex finding algorithms (see 
Section 5.2). The normalisation is to the number of entries and the simulation has been 
weighted to adjust the b-hadron production fractions, the B— >D branching ratios and the 
'wrong-sign' ^ Dg charm production rate to the same values as detailed in Section 8.1. 
The discrimination attained between signal and background is also plotted together with 
the distribution shape expected from tracks in u,d,s,c-events. The background distribu- 
tion in Figure 4 is dominated by tracks from weak b-hadron decays and the small spike at 
values close to BD-Net= —1 is due to scmi-leptonic decays of b-hadrons which arc read- 
ily recognised by the network as not coming from the D-vertex. The agreement between 
data and simulation is good and the effect of any residual discrepancy on the analysis 
was found to be insignificant (see Section 8.1). 




BD-Net 

Figure 4: Distributions of the BD-Net from 1994 data and simulation, normalised to the 
number of entries. 



Charm quarks produced from the upper- or W-vertex in b-hadron decay via e.g. W"'" — » cs. 
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5 Proper Decay Time Reconstruction 

This section deals with the reconstruction of the proper time defined as, 

t — Lmo/pc 

where L and p are the reconstructed decay length and momentum of the B-candidate 
respectively and mo is the B rest mass which was taken to be 5.2789 GeV/c^. The 
magnitude of the B-candidate momentum vector was fixed by the relationship E"^ = 
+ 771^, where E is the reconstructed B-candidate energy. The reconstruction of the two 
essential components of the proper time, namely the B-candidate energy and the decay 
length, are now dealt with in some detail. 

5.1 B-candidate Energy Reconstruction 

A novel method to reconstruct the B-candidate energy was used which involved train- 
ing a neural network (EB-Net) to return a complete probability density function (p.d.f.) 
for the energy on a hemisphere-by-hemisphere basis. Training the network proceeds by 
dividing up the inclusive truth B-energy distribution from the Monte Carlo f{B) into 
many equal slices or threshold levels. One neural network output node is assigned to 
each level and the target value for each training event becomes a vector with as many 
elements as there are output nodes, defined by the classification: 'is the energy of the B 
in this hemisphere, higher (1) or lower (-1) than the associated threshold at this output 
node'. The median of the resulting estimate of f{B) is taken as the EB-Net estimator 
for the B-energy. 

For the network training sixteen input variables were used which included different 
estimators of the energy available in the hemisphere together with some measures of the 
expected quality of such estimators e.g. as given by such quantities as hemisphere track 
multiplicity and hemisphere reconstructed energy. The key inputs to the network were: 

• In 2-jet events the sum (over all particles i in a hemisphere) of the vectors 
(pj, Ei^rrii) was formed, weighted by the TrackNet value (if particle i was a charged 
particle) or weighted by a function of the rapidity (if particle i was a neutral cluster). 
In this way particles from the b-hadron decay received a higher weight in the sum 
and hence an estimate of the b-hadron vector was obtained: Pb = (pb, -E_b, m^). 
For > 3-jet events Pb was estimated via the Rapidity Algorithm. These estimates 
of Eb were then used directly as an input to the EB-Net and the momentum vec- 
tors Pb provided the direction constraint for the secondary vertex fitting algorithms 
described in Section 5.2. 

• The estimate of the b-hadron energy Eb was corrected to account for sources of 
missing energy. The correction procedure was motivated by the observation in sim- 
ulation of a correlation between the energy residuals AE = Eb — E^"'^^"'^'^'^ and 
rriB, which is approximately linear in m^, and a further correlation between AE 
and Xh (the fraction of the beam energy in the hemisphere) resulting from neutral 
energy losses and inefficiencies. The correction was implemented by dividing the 
data into several samples according to the measured ratio Xh and for each of these 
classes the B-energy residual AE was formed as function of m^. The median values 
of AE in each bin of were calculated and their ms-dependence fitted by a third 
order polynomial 

AE{mB, Xh) = a-^ b{mB - {rriB)) + c(ms - {mB)f + d{mB - {mB)f 
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Figure 5: (a) The EB-Net and (b) the reconstructed B-candidate decay length residual 
i.e. the difference between the reconstructed value and the generated value based on 1994 
simulated data. 



The four parameters a, 6, c, d in each Xh class were then plotted as a function of Xh 
and their dependence fitted with second and third-order polynomials. Thus a smooth 
correction function was obtained from the simulation describing the dependence of 
AE on niB and on the hemisphere energy. The corrected b-hadron energy was the 
most powerful input variable to the EB-Net, having a correlation of 73% to the true 
b-hadron energy. 

The performance of the EB-Net estimator is shown in Figure 5(a) which plots the 
residual of the EB-Net variable with the generated energy value. A double-Gaussian fit 
to the distribution gives a central, narrow, Gaussian covering 67% of the total area with 
standard deviation of 2.5 GeV. (Note that the fits are approximate and are only to gauge 
the widths of the distributions. They are not used in the lifetime measurements.) 

5.2 Decay Length Reconstruction 

Starting from the standard secondary vertex described in Section 4.3, four algorithms 
were implemented, based on the BD-Net, with the aim of improving the decay length 
resolution and minimising any bias of the type described in Section 4.5, resulting from the 
inclusion of tracks from the cascade D-decay vertex in the B-decay vertex reconstruction. 
In addition to passing the standard quality cuts listed in Section 4.3, tracks were required 
to have TrackNet values > 0.5 to be considered for any of the four algorithms. 

1) In the Strip-Down method candidate tracks were selected if, in addition to the cuts 
described above, they had BD-Net values < 0.0. A secondary vertex fit was made 
if there were two or more tracks selected. If the fit failed to converge (under the 
same criteria as were applied to the standard fit - see Section 4.3), and more than 
two tracks were originally selected, the track with the highest contribution was 
removed and the fit repeated. This procedure continued iteratively until convergence 
was reached or less than two tracks were left. Technically, the fit was the same as 
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that used to fit the standard secondary vertex described in Section 4.3 except the 
starting point was the secondary vertex position estimate of the standard fit. 

2) In the D-Rejection method, a cascade D-candidate vertex was first built by fitting 
a common vertex to the two tracks with the largest BD-Net values in the hemisphere. 
If the invariant mass of the combination was below the D-meson mass, an attempt 
was made to include also the track with the next largest BD-Net value. This process 
continued iteratively until either the mass exceeded the D-mass, there were no further 
tracks, or the fit failed to converge. The B-candidate vertex was then fitted using 
the Strip-Down algorithm but applied to all tracks except those already selected for 
the D- vertex. 

3) In the Build-Up method those two tracks with TrackNet bigger than 0.5 and small- 
est BD-Net values were chosen to form a seed vertex. If the invariant mass of all 
remaining tracks with TrackNet > 0.5 exceeded the D-mass, that track with the 
lowest BD-Net output was also fitted to a common vertex with the two seed tracks. 
This process continued iteratively until either the fit failed to converge or the mass 
in remaining tracks dropped below the D-meson mass. 

4) The Semileptonic algorithm was designed to improve the vertex resolution for 
semileptonic decays of b-hadrons where energy has been carried away by the asso- 
ciated neutrino. When there was a clear lepton candidate in the hemisphere, the 
algorithm reconstructed a cascade D-candidate vertex in a similar way to the D- 
rejection method but with the lepton track excluded. The tracks associated with 
the vertex were then combined to form a 'D-candidate track' which was extrapo- 
lated back to make the B-candidate vertex with the lepton track if the opening angle 
between the lepton and D-candidate satisfied | cos0l£)| < 0.99. 

The choice of decay length for the decay time calculation was dictated by optimising 
the resolution and minimising any bias while still retaining the best possible efficiency. If 
more than one of the four algorithms was successful in reconstructing a vertex, the choice 
was made in the following order: 

1) the Strip-Down method was chosen if the algorithm had a decay length error smaller 

than 1 mm, 

2) if the Strip-Down method criteria were not met, the D-rejection method was used if 
the decay length error was smaller than 1 mm, 

3) if the criteria for 1) and 2) were not met the Build-Up vertex was chosen if the decay 
length error was smaller than 200 /im, 

4) if the criteria for 1), 2) and 3) were not met the Semileptonic algorithm was used if 
the decay length error was smaller than 1 mm. 

About one third of all hemispheres, passing the event selection cuts, were rejected by the 
decay length selection procedure in data and in simulation. There were 180010 vertices 
selected in the 1994 data set and 86796 in 1995. Figure 5(b) plots the residual between 
the reconstructed decay length and the generated value. A double-Gaussian fit to the 
distribution gives a central, narrow, Gaussian covering 71% of the total area with a 
standard deviation of 300 fim. The lack of a significant positive bias to this distribution, 
illustrates that the infiuence of tracks from cascade D-meson decays has been successfully 
minimised by employing algorithms based on the BD-Net. 
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6 Selection of B+ and B° Enhanced Samples 

The enrichment of B° and mesons was part of a general attempt to provide a 
probabihty for an event hemisphere to contain a b-hadron of a particular type. The 
result was implemented in a neural network {NN{Bx)) consisting of 16 input variables 
and a 4-node output layer. Each output node delivered a probability for the hypothesis 
it was trained on: the first supplied the probability for Bg mesons to be produced in the 
hemisphere, the second for B'^ mesons, the third for charged B-mesons and the fourth for 
all species of b-baryons. The method relied heavily on the reconstruction of the following 
quantities: 

• b-hadron type probabilities P(Bx): supplied by an auxiliary neural network con- 
structed to supply inputs to the more optimal NN{Bx) network. In common with 
the NN{Bx), there were four output nodes trained to return the probabihty that the 
decaying b-hadron state was B+, B°, B^ or b-baryon. There were fifteen input vari- 
ables in total, the most powerful of which were the hemisphere TrackNet-weighted 
charge sum, which discriminates charged from neutral states, and variables that 
exploit the presence of particular particles produced in association with b-hadron 
states. Examples of this include B^ mesons, which are normally produced with a 
charged kaon as the leading fragmentation particle with a further kaon emerging 
from the weak decay, and in B+ and B° production where the decay is associated 
with a larger multiplicity of charged pions than that for B^ and b-baryons (which 
on average will produce a higher proportion of neutrons, protons and kaons). 

• The b-hadron flavour i.e. the charge of the constituent b-quark both at the frag- 
mentation (Ffrag.) and decay time (Fdec.)- knowledge of b-hadron flavour provides 
the network with valuable information about whether a B° state was present since 
the fragmentation and decay flavour will, on average, disagree for the case where the 
B° oscillated. The flavour was determined by flrst constructing, with neural network 
techniques, the conditional probability for each charged particle in the hemisphere to 
have the same charge as the b-quark in the b-hadron. This was repeated separately 
for each of the four possible b-hadron type scenarios i.e. B°,B"'",Bg or b-baryon. The 
flavour network was trained on a target value of 1) if the particle charge was 
correlated(anti-correlated) to the b-quark charge. The main input variables were 
those related to the identification of kaons, protons, electrons and muons together 
with quantities sensitive to the B-D vertex separation in the hemisphere. Tracks 
originating from the fragmentation (decay) phase were discriminated by checking 
that the TrackNet value is less (greater) than 0.5. In a final step, these track level 
probabilities were combined via a likelihood ratio into hemisphere quantities. Pro- 
viding separate flavour networks for the different b-hadron types not only ensured 
that the information was optimal for the case of B° but also helped the performance 
of the enrichment NN{Bx) network by providing information that was speciflc to a 
particular b-hadron type. 

The inputs to the NN{Bx) were constructed to exploit optimally all of the information 

that the b-hadron production and decay process reveals. The basic construct for input 
variables was the following combination of the flavour and b-hadron type information 
described above: 

Fdeci^x) ■ Ffrag.(Bx) ' Pi^x)- (1) 

The upper plots of Figure 6 show the output of the B"*" and B° output nodes of the 
NN[Bx) in simulation and data. The simulation is further divided into the different 
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Figure 6: The upper plots show the output of the and B° output nodes of the 
NN{Bx) in the 1994 data and simulation for the different b-hadron types. The nor- 
malisation is to the number of data events and overlaid is the b-hadron composition as 
seen in the Monte Carlo. The lower plots trace the change in purity of the different 
b-hadron types, per bin, as a function of cuts on the NN{B~^) and NN{B^) respectively. 



b-hadron types and the lower plots trace the change in purity of the different types in 
each bin of the network output at the B"*" and B*^ output nodes respectively. Note that 
the background from u,d,s- and c-events is labelled as 'bg'. 



7 Extraction of and B Lifetimes 

This section describes how the data and simulation samples were prepared, gives de- 
tails of the fitting procedure itself and summarises the results obtained. Section 7.1 hsts 
corrections made to the default simulation to account for known discrepancies with data 
and to update b-physics parameters to agree with recent world measurements. In Sec- 
tion 7.2 the final selection and composition of the B"*", B° and Tb samples are described 
together with an explanation of how the region at low proper times i.e. < 1 ps was 
handled in the hfetime fits. Section 7.3 gives technical details of the fit and presents 
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the results obtained. Lifetimes were measured separately in 1994 and 1995 data and then 
combined to give the final results which are presented in Section 9. 

7.1 Simulation Weighting 

Weights were applied to the simulation to correct for the following effects: 

• The world average of measurements of B^, Ab lifetimes and b-hadron production 
fractions, as compiled for the Winter conferences of 2001 [13], listed in Table 1. 
Note that using more [13] recent world average values has a negligible effect on the 
analysis results. 

• The Peterson function used in the Monte Carlo {{x) = 0.706) was changed to agree 
with the functional form unfolded from the 1994 DELPHI data set [12] {{x) ~ 
0.7153). 

• A hemisphere 'quality flag' which was proportional to the number of tracks in the 
hemisphere likely to be badly reconstructed e.g. those tracks failing the stan- 
dard selection criteria of Section 4.3. The weight was constructed to account for 
data/simulation discrepancies in this variable and was formed in bins of the number 
of charged particles in the hemisphere, that passed the standard quality cuts of Sec- 
tion 4.3, thus ensuring that, overall, the multiplicity of good charged particles was 
unchanged after the application of the weight. 



b-hadron species 


Lifetime 


Production fractions /(& — > Bx) 


Values 




p{Ba:, baryon) 




1.464 ± 0.057ps 


/(B^) =0.097 ±0.011 




+0.034 


b-baryons 


1.208 ± O.OSlps 


/(baryon) =0.104 ±0.017 


+0.034 




or B+ 




/(Bd) = /(B+) = 0.399 ± 0.010 


-0.577 


-0.836 



Table 1: Values for the b-hadron lifetimes and production fractions (together with corre- 
lations) used to re-weight the Monte Carlo. 



7.2 Fit working point 

The selection conditions imposed on the data samples used for the lifetime fits were 
motivated by the wish to minimise the total error on the final results. Systematic error 
contributions due to inexact detector resolution simulation and the physics modelling of 
u,d,s and charm production, imply that relatively high b-hadron purities were required 
while still keeping the selection efficiency above a level where the statistical error would 
begin to degrade significantly. 

With these considerations in mind, the final data samples to be used in the fitting 
procedure were selected by cutting on the NN{Bx) neural network outputs, described 
in Section 6, at > 0.52 and > 0.6 respectively to obtain enhanced samples in B+ and 
B°. These cut values corresponded to a purity in both B"^ and B° of approximately 70% 
according to the simulation. No reconstructed hemisphere passed both the B+ and B*^ en- 
hancement cuts simultaneously and hence the two samples were statistically independent. 

The region below about 1.0 ps in proper lifetime is particularly challenging to simu- 
late. The modeUing of very small hfetimes is rather sensitive to details of reconstruction 
resolution and the modelling of events which contain no intrinsic lifetime information 
such as u,d,s events and the reconstruction or spurious vertices. In addition, the lower 
plots of Figures 7 and 8 show that the purities of the different b-hadron types is rapidly 
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changing in this region, making them particularly difficult to model. These issues meant 
that the low lifetime region was not well enough under control systematically for precision 
lifetime information to be extracted and the region below 1.0 ps was therefore excluded 
from the analysis. This point is illustrated in Figure 10 which shows that the fit results 
only become stable in all samples for a fit starting point larger than 1.0 ps. 

After all selection cuts already described, the size and composition of the B+and 
B° enhanced samples are summarised in Table 2. The B"'"(B°) sample sizes correspond 
to a selection efficiency, with respect to the starting number of B"*" and B*^ states in the 
hadronic sample, of 10.1% (3.8%) for both 1994 and 1995 data. 





B+ Sample 


B" Sample 


Real data sample size 1994(1995) 


27356(13150) 


9293(4335) 


Simulation sample size 1994(1995) qq 


61821(23415) 


22667(8533) 


Simulation sample size 1994(1995) hh 


161198(42951) 


59015(15863) 


B+ fraction 


71.7% 


15.8% 


B" fraction 


20.7% 


68.5% 


Eg fraction 


3.7% 


10.1% 


b-baryon fraction 


2.4% 


4.4% 


u,d,s fraction 


0.2% 


0.1% 


c fraction 


1.2% 


1.0% 



Table 2: The B"*" and B° sample size in data and simulation and the composition of the 
simulation. The simulation has been weighted for the quantities listed in Section 7.1. 

The data sample used to fit for the mean b-hadron lifetime passed through the same 
event selection and proper time cuts as for the B+ and B° samples but without any 
requirement on the NN{Bx) neural network outputs. The size and composition of the 
mean b-hadron sample is summarised in Tabic 3. These numbers imply that the mean 
lifetime measurement is valid for a b-hadron mixture, as given by the simulation, of 
B+= 41.9%, B0= 41.2%, B0= 8.9% and b-baryon= 8.0%. 





b-hadron 


Sample size 1994(1995) 


114317(54958) 


Simulation sample size 1994(1995) qq 


262697(98230) 


Simulation sample size 1994(1995) bb 


677998(180634) 


B+ fraction 


41.2% 


B° fraction 


40.5% 


Bg fraction 


8.7% 


b-baryon fraction 


7.7% 


u,d,s fraction 


0.3% 


c fraction 


1.7% 



Table 3: The mean b-hadron sample size in data and simulation and the composition of 
the simulation. The simulation has been weighted for the quantities listed in Section 7.1. 



7.3 Lifetime Results 

The B"*", B° and Tb lifetimes were extracted by fitting the simulated proper time 
distribution to the same distribution formed in the data using a binned method. As 
discussed in Section 7.2 the start point of the fit range was chosen to be 1 ps. The upper 
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limit was positioned to avoid the worst effects of spurious, mainly two-track vertices, 
with very long reconstructed lifetimes while still accepting the vast majority of the data 
available. Nominally 100 bins were chosen but the exact binning was determined by the 
requirement that at least 10 entries be present in all bins of the data distribution. 

To avoid the need to generate many separate Monte Carlo samples with different B- 
lifetimes, weighting factors were formed for each lifetime measurement from the ratio of 
exponential decay probability functions. Specifically, the weight, 

Told f tii^Tnew '^old) 

-exp 



for measurement i and true B-lifetime tj, effectively transforms the Monte Carlo lifetimes 
generated with a mean lifetime Toid to be distributed with a new mean value of Tnew 
Throughout the fit for the B"*" and B° lifetimes, the Bg and lifetime components 
were weighted to the current world average numbers listed in Table 1 and for the Tb fit, 
the starting value in the simulation was 1.6 ps. The function given below was then 
minimised with respect to the B"*" and B° lifetimes in a simultaneous two parameter fit 
or to the mean b-hadron lifetime Tb in a one parameter fit. 



BO,B+ 



"bins ('\jTS'im '\Tdata\2 



Here, Nf"-*"- is the number of data entries in bin i and W?'^'^ is the corresponding sum of 
weights in bin i of the simulation. 

The results from all lifetime fits, after imposing the working point conditions and 
following the above procedure, are listed in Table 4. In the table, the first error quoted 
is statistical and the second systematic. The various sources of systematic error are 
described in Section 8. Results are given for 1994 and 1995 data separately and combined 
taking into account correlated systematic errors as described in Section 9. 



b-State 


Fitted Lifetime 




'94 


'95 


Combined 


B+ 


1.624 ±0.017 ±0.023 ps 


1.623 ±0.025 ±0.019 ps 


1.624 ±0.014 ±0.018 ps 


B» 


1.548 ±0.026 ±0.035 ps 


1.497 ±0.039 ±0.041 ps 


1.531 ±0.021 ±0.031 ps 


'b+ 

Tp,0 


1.049 ±0.025 ±0.027 


1.085 ±0.040 ±0.036 


1.060 ±0.021 ±0.024 


Tb 


1.577± 0.006 ±0.008 ps 


1.555 ±0.009 ±0.011 ps 


1.570 ±0.005 ±0.008 ps 



Table 4: The results of the lifetime fits in the 1994 ^995 data samples where the first 
error quoted is statistical and the second systematic. 



The B"*" and B° fits (to the 1994 data) are shown in Figure 7. The correlation coefficient 
between B"*" and B'^ hfetimes was found to be —0.51 for both the 1994 and 1995 fits. The 
fit at the minimum point was 178 for 160 degrees of freedom for 1994 data and 142 for 
143 degrees of freedom for 1995 data. 

The mean b-hadron lifetime fit is shown in Figure 8. The the minimum point 

was 76 for 88 degrees of freedom in 1994 data and 70 for 88 degrees of freedom for 1995 
data. 

Figure 9 illustrates the effect of cut scans in the Z° bb event purity (i.e. the 
estimated fraction of the data sample fitted coming from Z° bb events) showing a 
good stability over a wide range of the cut values for the lifetime ratio t{B^)/t{B^) 
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Figure 7: The upper two plots show the result of the fit in the B'^ (left) and (right) 
samples in 1994 (histogram) compared to data(points). The b-hadron composition of the 
andB^ sample is also indicated where 'bg' refers to the background from non-hh TP de- 
cays. The lower two plots trace how the fractional composition of the sample changes in 
bins of the reconstructed lifetime. The vertical line at t = 1 ps indicates that data below 
this point are removed from the analysis. 



and the mean b-hadron hfetime. Similarly Figure 10 illustrates that the results are very 
stable over a wide range of different start points for the fits above the default cut point 
of t = 1 ps. 

A further crosscheck on the results was made by repeating the fits for one data set 
using the simulation sample compatible with another data set e.g. fitting 1994 data using 
1995 simulation. It was found that all fit results (for B+, B° and Tb) for both cases (1994 
data using 1995 simulation and 1995 data using 1994 simulation) changed by amounts 
that were within the systematic error for detector effects quoted in Table 5 which provides 
a rough check that aspects of detector and physics modelling are well under control. 
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Figure 8: The left plot shows the result of the mean h-hadron lifetime fit in 1994 (his- 
togram) compared to the data(points) at the working point. The right plot traces how the 
fractional composition varies in bins of the reconstructed lifetime. The vertical line at 
t — l^s indicates that data below this point are removed from the analysis. 



8 Systematic Uncertainties 

Systematic uncertainties on the lifetime measurements come from three main sources. 
The first source is from the modelhng of heavy flavour physics parameters in our Monte 
Carlo generator. Since attempts were made to model these effects using current world 
averages, these errors are largely irreducible. The second source comes from the analysis 
method itself and the choices made in determining the measurement working point. The 
good level of agreement between simulation and data and the fact that the result is 
stable within a wide range of the working point (e.g. as shown in Figure 9) mean that 
these errors are kept to a minimum. The third source of systematic uncertainty can 
be generically termed 'detector effects' and results from a less-than-perfect modelling in 
simulation of the response of the detector. 

Tables 5 and 6 present the full systematic error breakdown for the measurements of 
r(B+), r(BO) and in 1994 and 1995 data. 

8.1 Heavy Flavour Physics Modelling 

Where possible, B-physics modelling uncertainties were estimated by varying central 
values by plus and minus one standard deviation and taking half of the observed change 
in the fltted lifetime value as the resulting systematic uncertainty from that source. 

The Bg, b-baryon lifetimes, b-hadron production fractions and fragmentation (x) value 
have been varied within their errors as listed in Section 7.1 and half of the full variation 
in the results has been assigned as an error. In the case of the b-hadron production 
fractions, the variation was made taking into account correlations from the covariance 
matrix listed in Table 1. 

Close attention was paid to possible systematic effects on the analysis due to the mod- 
elling of b ^ charm branching ratios, where current experimental knowledge is scarce. 
The charm content impacts on the performance of the B+ and B° enhancement networks 
and can pull the reconstructed B- vertex position to longer decay lengths. The size of this 





rB+ (1994|1995) 


tbo (1994|1995) 


^ (1994|1995) 


Result [ps] 


1.6241 


1.6233 


1.5483 


1.4971 


1.0492 


1.0848 


Statistical Error [ps] 


0.0168 


0.0251 


0.0255 


0.0388 


0.0247 


0.0398 


Source of Systematic Error 


Range 


Atb+ [ps] 


Atbo [ps] 




Physics Mock^Uiug 










Bg lifetime 


f .464 ± 0.057 ps 


0.0007 


0.0080 


0.0050 


b-baryon lifetime 


1.208 ± 0.051 ps 


0.0007 


0.0030 


0.0028 


b-hadron prod, fractions 


O J. J. 

bee text 


0.0035 


0.0035 


0.0004 


Fragmentation function 


\X) = U. ( ioo ± U.UUOZ [izj 


0.0037 


0.0026 


0.0040 


B DX branching fractions 


See text 


0.0081 


0.0086 


0.0017 


BR(B — > wrong-sign charm) 


20.0% ± 3.3% 


0.0030 


0.0047 


0.0014 


BR(Bs ^ D,) 


35% ^ 70% 


0.0017 


0.0076 


0.0062 


BR(b ^c-baryon X) 


9.6% ± 3.0% 


0.0009 


0.0032 


0.0016 


D+,D° topo. branching ratios 


[14] 


0.0016 


0.0113 


0.0084 


B meson mass 


ms = 5.2789 ± 0.0018 Gev/c^ 


0.0004 


0.0015 


0.0011 


b-hadron Reconstruction 










b and c efficiency correction 


On/off 


0.0076 


0.0048 


0.0050 


0.0066 


0.0015 


0.0017 


NN{B^) cuts 


65%-75% purity 


0.0093 


0.0133 


0.0216 


0.0255 


0.0196 


0.0269 


NN{B^) shape 


See text 


0.0008 


0.0006 


0.0099 


0.0126 


0.0059 


0.0097 


Sec. vertex multiplicity 


See text 


0.0022 


0.0019 


0.0042 


0.0040 


0.0040 


0.0042 


Detector Effects 










Resolution and 


On/off 


0.0171 


0.0075 


0.0144 


0.0200 


0.0110 


0.0159 


hemisphere quality 
















Total Systematic Error 




0.0234 


0.0192 


0.0345 


0.0406 


0.0269 


0.0355 
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Figure 9: The variation in the fitted lifetimes as a function of the bb purity for the 

ratio t{B^)/t{B^) in (a) 1994 ^^'^ (i') 1995 data and for the mean b-hadron lifetime in 
(c) 1994 and (d) 1995 data. The upper and lower shaded hands represent the statistical 
one standard deviation errors which are correlated bin-to-bin. 



pull in turn depends on whether a D'' or D+ was produced since t{D~^) ~ 2.5 times larger 
than t(D°). Specific aspects of the Monte Carlo that were found to warrant systematic 
error contributions were: 

(a) Inclusive, B — » D branching ratios were adjusted in the Monte Carlo according to a fit 
using all currently available measurements from [2] as constraints. The values taken were: 
BR {B° ^D+X) = 15.6%, BR{B'^ ^B^X) = 65.8%, (B" D+X) = 29.3%, 
BR{B^ D°X) = 52.1%. The full difference seen in the results due to this change 
was assigned as a systematic error. 

(b) The standard Monte Carlo data set used contained a wrong-sign charm production 
rate into Dg of 11%. This rate is now known to be too low due to the production of 
D^/D~^ mesons at the W vertex (in addition to Dg), and an estimate of the overall rate 
based on [15] is BR(wrong-sign)= 20.0%±3.3%. To account for this discrepancy with cur- 
rent measurements, the wrong-sign rate in the simulation was weighted up to 20.0% and 
all quoted lifetime results were shifted to be valid for this higher wrong-sign rate. An 
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Figure 10: Lifetime fit results as a function of varying the start point of the fit for B"*" in 
(a) 1994 and (b) 1995, B° in (c) 1994 (d) 1995 and for the mean b-hadron fit in (e) 
1994 and (f) 1995. 
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TB for 1994 


TB for 1995 


Result [ps] 


1.5773 


1.5546 


Statistical Error [ps] 


0.0059 


0.0085 


Source of Systematic Error 


Range 


Atb [ps] 


Atb [ps] 


Physics Modelling 








b-hadron prod, fractions 
Fragmentation function 
BR (B wrong-sign charm) 
BR(Bs ^ Ds) 
BR(b ^c-baryon X) 
B — > DX branching fractions 
D+ , D" topo. branching ratios 
B-meson mass 


See text 
(a;) = 0.7153 ± 0.0052 [12] 
20% ± 3.3% 
35% 70% 
9.6% ± 3.0% 
Sec text 
See text 
ruB = 5.2789 ± 0.0018 Gev/c^ 


0.0011 
0.0038 
0.0010 
0.0012 
0.0017 
0.0008 
0.0008 
0.0004 


b-hadron Reconstruction 








b and c efficiency correction 
Sec. vertex multiplicity 


On/off 
See text 


0.0057 
0.0019 


0.0048 
0.0023 


Detector Effects 








Resolution and 
hemisphere quality 


See text 


0.0035 


0.0083 


Total Systematic Error 




0.0084 


0.0109 



Table 6: Summary of systematic uncertainties in the mean b-hadron lifetime for 1994 
and 1995 data. Systematic errors are assumed independent and added in quadrature to 
give the total systematic error. 



error was then assigned based on the change in the hfetime results observed when the 
wrong-sign rate was further changed by ±3.3%. The impact of the simulation contain- 
ing only wrong-sign Dg mesons instead of a mixture of Dg, D*^ and D+ was tested with 
specially generated Monte Carlo data sets, and the effect found to be small compared to 
the overall effect of having almost double the rate of events containing two D-mesons per 
hemisphere. 

(c) BR(Bs DgX) is currently known to, at best, ±30% [2] and was varied in the Monte 
Carlo by a factor two from the default value of 35%. The full change in the fitted lifetime 
was then assigned as a systematic error. 

(d) BR(b — >c-baryon X), where b represents the natural mixture of b- mesons and baryons 
at LEP, was varied from the default value in the simulation of 9.6% by ±3.0%. This 
range covers the uncertainty on this quantity from experiment which currently stands at 
BR(b — >c-baryon X)= 9.7 ± 2.9% [2]. Half of the full change in the fitted lifetime was 
then assigned as a systematic error. 

The uncertainty from D-topological branching fractions was estimated from the dif- 
ference in the fit result obtained when weighting according to the results from [14]. The 
masses of B-mesons were varied within plus and minus one standard deviation of the 
value assumed in the BSAURUS package and half of the change seen taken as a system- 
atic error. 

Since many of the physics modelling systematics investigated are significantly smaller 
than the statistical precision, the approach was taken to average the errors, evaluated in 
1994 and 1995 data, weighted by the statistical error for each year. This ensures that the 
effects of statistical fluctuations in the determination of these errors are minimised and 
explains why the physics modelling errors appearing in Table 5 are the same for 1994 and 
1995. 
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8.2 b-hadron Reconstruction 

The efficiency for reconstructing bb and cc events (as a function of the event b-tag) 
has been extracted from the real data by a double-hemisphere tagging technique. At 
the b-tag value of the working point, the results of this study suggest that while the 
reconstruction efficiency for bb events might be underestimated in the simulation by 
about ~ 2% (relative), the efficiency for cc events in simulation is ~ 8% (relative) lower 
than in data. To account for this possible source of error the difference seen in the fit 
results, when these efficiencies were changed in the simulation to agree with the numbers 
above, was assigned as a systematic error. Since a large part of the discrepancy between 
simulation and data in the cc event reconstruction efficiency is due to a less-than-perfcct 
modelling of charm physics, this error contribution has already been partially accounted 
for by the explicit charm physics systematics detailed above. Given the current level of 
uncertainty in this sector, the conservative approach of quoting both error contributions 
is preferred. 

As was remarked in Section 7.3, uncertainties resulting from the method itself have 
been checked by scanning regions around critical cut values to check for stability as 
illustrated in Figure 9. In addition, the binning used for the formulation was varied 
over a wide range as was the minimum number of entries per bin (set by default to 10) 
and were both found to give no significant change in the results. 

The impact on the analysis of any residual discrepancy between data and simulation 
in the BD-Net variable was checked by studying the effect of removing the Strip-Down 
vertex algorithm (see Section 5.2) from the analysis. Of the four vertex algorithms used, 
the Strip-Down method is the most sensitive to details of the BD-Net variable since it 
imposes a direct cut in the BD-Net as part of the track selection. Removing the Strip- 
Down algorithm and replacing it with one of the other three methods, selected by the 
same criteria as described in Section 5.2, resulted in lifetime results for 1994 data that 
changed by At{B+) = -0.0017 ps and At(S°) = -0.0063 ps i.e. weU within the total 
systematic error quoted. In addition it was confirmed that the proper time distributions 
were well compatible when the cut imposed on the BD-Net distribution was changed from 
the default value of zero to ±0.2. 

The analysis assumes that the B+ and B° purities are well modelled by the simulation. 
A systematic error will arise if the shape of the NN{B^) and NN[B^) network outputs 
differ from the data and/or the composition of the B+ and B*^ simulated samples differ. 
Any difference in shape between data and simulation was accounted for in the following 
way. It was assumed that the difference could be wholly accounted for by a change in just 
the B+ composition for the case of the NN{B'^) network and by just the B° composition 
for the case of the NN{B^) network. In this way it was found that the maximum error 
made by assuming that the B"^ and B° purities in the simulation were correct, was of 
order 2% and 4%, respectively. The effect of these changes were then propagated into 
errors on the extracted lifetimes. To account for any composition differences, half of the 
maximum variation in the fitted lifetime while scanning the purity range [65%, 75%] was 
assigned as an error. The scan range was chosen to enclose the largest uncertainty on the 
purity found from analysing the shapes of the network outputs described above. 

The multiplicity of tracks in the reconstructed b-hadron vertex was found to be in 
overall good agreement between the data and simulation. To account for any residual 
differences a weight was formed from the ratio of the data and simulation distributions 
and the change seen as a result of applying this weight was assigned as a systematic error. 



24 



8.3 Detector Modelling 

In order to account for uncertainties in the simulation originating from detector re- 
sponse modelling, the effect of switching on and off the following corrections was studied: 

• the hemisphere quality weight, described in Section 7.1, 

• an attempt to obtain a better match of the track impact parameter and error (with 
respect to the primary vertex) between simulation and data according to the pre- 
scription detailed in [7]. 

Since in general, knowledge of detector modeUing uncertainties are not at the same level of 
understanding as e.g. our knowledge of the difference between the B-production fractions 

in our Monte Carlo and the world averages, the following approach was taken to assigning 
systematic values for these effects: all four combinations of switching these corrections 
on/off for the B"'"/B° fit and the Tb analysis were made and the fitted lifetimes of the four 
possibihties recorded. The central results for B+,B° and Tb were then chosen to be the 
mean values of these four combinations, and the resulting systematic error from detector 
response modelling was assigned to be half of the maximum spread of the values recorded 
from the four combinations. This error is listed in tables 5 and 6 as 'Resolution and 
hemisphere quality'. 

8.4 Closing Remarks on Systematic Errors 

In general it can be concluded from Tables 5 and 6 that detector effects dominate. 
Physics modelling errors come essentially only from b-physics sources, since the contami- 
nation from light-quark and charm events is so small, and are generally well under control. 
For the case of the B+ and B'^, additional systematic error contributions arising from the 
enhancement neural networks ( NN{Bx)) reflect the difficult task of modelling accurately 
such complex variables. 

9 Summary and Conclusion 

The lifetimes of B+, B°, their ratio and the mean b-hadron lifetime have been mea- 
sured. The analysis isolated b-hadron candidates with neural network techniques trained 
to exploit the physical properties of inclusive b-hadron decays. Binned fits to the 
resulting DELPHI data samples collected in 1994 and 1995 yielded the results presented 
in Table 5 for B+ and B° and the result for the mean b-hadron lifetime is presented in Ta- 
ble 6. The results for 1994 and 1995 were combined, treating all systematic contributions 
as independent apart from the following (which were assumed to be 100% correlated): 

• all physics modelling errors, 

• the NN{Bx) shape error, 

• the secondary vertex multiplicity error. 

The combined results for the B"*" and B° were, 

tb+ = 1.624±0.014 (stat)±0.018 (syst) ps 
tbo = 1.531±0.021 (stat)±0.031 (syst) ps 
^ = 1.060±0.021 (stat) ±0.024 (syst) 

and for the average b-hadron lifetime: 

Tb = 1.570±0.005 (stat)±0.008 (syst) ps. 
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(Note that the average b-hadron hfetime result is vahd for a b-hadron mixture, as given 
by the simulation, of B+= 41.9%, B°= 41.2%, 3°= 8.9% and b-baryon= 8.0%). 

These results are well compatible with previous DELPHI results in this area using the 
1991-1993 data sets: based on D^*H reconstruction [16] and [17] and from an inclusive 
secondary vertex approach [18]. No attempt has been made to combine these older results 
with the current analysis because of the vast difference in precision e.g. the error on the 
hfetime ratio ^ is now a factor five better than was achieved in [18]. 

Compared to existing measurements, the B+ lifetime result is currently the most accu- 
rate and is well compatible with all other measurements and with the world average value 
of 1.671 ±0.018 ps [2]. The precision of the B*^ lifetime result is similar to the best so far 
achieved from Z° data i.e. from an OPAL analysis based on inclusive D'^*H reconstruc- 
tion (1.541 ± 0.028 ± 0.023 ps [19]) and to recent results from the B-factory experiments, 
BABAR [20] and BELLE [21]. In addition, the B° hfetime result is well compatible with 
all other measurements and with the current world average value of 1.537 ±0.015 ps [2]. 
All published measurements of the lifetime ratio are presented in Table 7 together 

with their average [2]. 



Experiment 


Method 


Data set 


Ratio tb+ I tbo 


Reference 


ALEPH 




91-95 


1.085±0.059± 0.018 


[22] 


ALEPH 


Exclusive rec. 


91-94 


-1 97+0.23+0.03 
-0.19-0.02 


[23] 


CDF 




92-95 


1.110±0.056+°°^i^ 


[24] 


CDF 


Excl. (J/i'K) 


92-95 


1.093 ± 0.066 ± 0.028 


[25] 


DELPHI 


dw^ 


91-93 


1.00l[5-i^± 0.10 


[16] 


DELPHI 


Charge sec. vtx. 


91-93 


1.06l°'j? ± 0.10 


[18] 


L3 


Charge sec. vtx. 


94-95 


1.09± 0.07 ± 0.03 


[26] 


OPAL 




91-93 


0.99±0.14l°-°^ 


[27] 


OPAL 


Charge sec. vtx. 


93-95 


1.079± 0.064 ± 0.041 


[28] 


SLD 


Charge sec. vtx. £ 


93-95 


1.03+°:i^± 0.09 


[29] 


SLD 


Charge sec. vtx. 


93-95 


1.0ll[5-[;^± 0.05 


[29] 


BABAR 


Exclusive rec. 


99-00 


1.082 ± 0'.026± 0.012 


[30] 


BELLE 


Exclusive rec. 


99-01 


1.091 ± 0.023± 0.014 


[21] 


Average 






1.085 ± 0.017 





Table 7: Measurements of the ratio tb+/tbo. 



It can be seen that the result from this analysis is currently the most precise of the 
measurements from Z° decay data and the CDF experiment at the Tevatron and also 
has a precision comparable to the B-factory experiments BABAR and BELLE. Within 
the quoted errors the result is also compatible with all measurements and with the world 
average value. 

The result for the mean b-hadron lifetime significantly improves on the most precise 
existing measurement from L3, Tb = 1.556±0.010±0.017 ps [32] and is in good agreement 
with the most precise previous DELPHI publication on this subject [33]. In addition it is 
compatible with the current world average, 1.564 ±0.014 ps [2], which has been compiled 
assuming that all measurements are based on b-hadron samples with the same mixture 
of b-hadron species i.e. the b-hadron production fractions from decay. It is also infor- 
mative to compare these values with the inclusive lifetime defined as. 
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Tb = /(Bi)r(Bj), calculated using the current world average values for b-hadron pro- 
duction fractions /(Bj) and lifetimes T(Bj) from [2]: Tb = 1.543 ps. 
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